Analysis of Parametric Oscillatory Instability 
in Signal Recycled LIGO Interferometer 
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We present the analysis of undesirable effect of parametric oscillatory instability in signal recycled 
LIGO interferometer. The basis for this effect is the excitation of the additional (Stokes) optical 
mode, with frequency u>i , and the mirror elastic mode, with frequency u> m , when optical energy 
stored in the main FP cavity mode, with frequency wo, exceeds the certain threshold and the 
frequencies are related as tuo — on + cu m . We show that possibility of parametric instability in this 
interferometer is relatively small due to stronger sensitivity to detuning. We propose to "scan" the 
frequency range where parametric instability may take place varying the position of signal recycling 



I. INTRODUCTION 

The full scale terrestrial interferometric gravitational 
wave antennae LIGO arc working now and have sensi- 
tivity, expressed in terms of the metric perturbation am- 
plitude, only 2 or 3 times worse than the planned level 
of h. ~ 1 x 10~ 21 0,0 (see current sensitivity curve in 
0|) • In Advanced LIGO (to be realized in approximately 
2012), after the improvement of the isolation from noises 
in test masses (the mirrors of the 4 km long optical FP 
cavities) and increasing the optical power circulating in 
the resonator up to W ~ 830 kW the sensitivity is ex- 
pected to reach the value of h. ~ 1 x 1 0~ 22 0, [B| . 

In we have analyzed the undesirable effect of para- 
metric oscillatory instability in the Fabry-Perot (FP) cav- 
ity, which may cause a substantial decrease of the anten- 
nae sensitivity or even the antenna malfunction. This 
effect appears above the certain threshold of the optical 
power W c circulating in the main mode, when the dif- 
ference u>o — between the frequency cuq of the main 
optical mode and the frequency u>i of the idle (Stokes) 
mode is close to the frequency cu m of the mirror me- 
chanical degree of freedom. The coupling between these 
three modes arises due to the ponderomotive pressure of 
light in the main and Stokes modes and the paramet- 
ric action of mechanical oscillation on the optical modes. 
Above the critical value of light power W c the ampli- 
tude of mechanical oscillation rises exponentially as well 
as optical power in the idle (Stokes) optical mode. How- 
ever, E. D'Ambrosio and W. Kells have shown Q that if 
in the same unidimensional model the anti-Stokes mode 
(with frequency cui Q = cuo + cu m ) is taken into account, 
then the effect of parametric instability will be substan- 
tially dumped or even excluded. In , we have presented 
the analysis based on the model of power recycled LIGO 
interferometer and demonstrated that anti-Stokes mode 
could not completely suppress the effect of parametric 
oscillatory instability. As possible "cure" to avoid the 
parametric instability we have proposed to change the 
mirror shape and introduce low noise damping [2j. D. 
Blair with collaborators proposed valuable idea to heat 



the test masses in order to vary curvature radii of mir- 
rors in interferometer and hence to control detuning and 
decrease overlapping factor between optical and acoustic 
modes 0, ^2 113 • Recently, the parametric instability 
effect was observed in experiment [lCj . 

It is interesting that the effect of parametric instability 
is important not only for large scale LIGO interferometer 
but it was observed also for micro scale whispering gallery 
optical resonators 0, ^| . 

In this paper we propose the detail analysis of para- 
metric instability in signal recycled Advanced LIGO in- 
terferometer (i.e. with additional signal recycling (SR) 
mirror) and show that, on the one hand, the paramet- 
ric instability in this interferometer can appear at low 
optical power (about several Watts) but, on the other 
hand, the probability that parametric instability condi- 
tion will be fulfilled is small due to small relaxation rates 
of optical modes (about several Hz). We also show that 
by varying the position of SR mirror one can change the 
frequency of anti-symmetric optical mode of interferom- 
eter and hence to scan the range of frequencies where 
parametric instability or its precursors may arise. 

In scction[n]wc derive the parametric instability condi- 
tions in signal recycled LIGO interferometer, we discuss 
these results in section Hill The details of calculations we 
present in Appendices. 



II. SIGNAL RECYCLED INTERFEROMETER 

We consider the LIGO interferometer with signal re- 
cycling (SR) and power recycling (PR) mirrors — see 
fig. ^ and notations over there. The port Eg is used for 
signal detection. Interferometer is pumped via port F5. 
Simplifications are the following: 

• Optical losses in all mirrors are absent (notes about 
generalization for non zero losses are made in the 
end of Scc. lIIIf) . Suspension noise in mirrors is also 
absent. 

• Transmittivity of input mirrors and length of both 
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FIG. 1: Signal and power recycled LIGO interferometer. Here 
Fi , Ei are the field amplitudes on the input mirror of the FP 
cavity 1 (and on plane (i)), F2, E2 are the field amplitudes on 
the input mirror of the FP cavity 2 and on the beam splitter, 
F3, E3, F4, E4 are the field amplitudes on the beam splitter. 



FP cavities are the same. They are tuned in reso- 
nance with the main mode. 



• Optical power W circulating inside the arms of in- 
terferometer is constant (approximation of constant 
field). 



• The distances between the input FP mirrors and 
beam splitter, and between the beam splitter and 
PR, SR mirrors are short (about several meters) — 
hence we consider the phase advance of waves trav- 
eling between these mirrors as constant and omit 
its dependence on frequency. 



FP cavities in arms 



We denote the mean amplitude of the wave in the main 
mode (with frequency Wo) by cursive capital letters and 
small fluctuation amplitude of the wave in the Stokes 
mode (u)i) by small letters. For example, the electrical 
field of wave falling on back mirror in FP cavity 1 is the 



following: 
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Here W is the mean power in the main mode circulating 
inside the cavity, c is the speed of light, dimcnsionless 
functions *4oin (?_!_), AMn[rj_) describe the distributions 
of optical fields over the beam cross section, integration 
J dfj_ is taken over the mirror surface. For simplicity 
below we consider the distributions of optical fields to be 
identical at all four mirrors in the arms. 

We write down the displacement vector of the elastic 
mode with eigenfrequency tu m , for example, for the end 
mirror of the FP cavity 1 as a product of time and space 
dependent functions: 

(xi (t) e - itUmt + xf {t)e iw ^) u(f) 

where xi and x* are slowly varying amplitudes and u 
is spatial vector of displacements of elastic mode in the 
mirror. 

Here we assume that the input and the end mirrors 
in the FP cavity 1 are elastically identical (the modes 
of their elastic oscillations coincide). Then we start from 
the formulas derived in Appcndix^in frequency domain: 
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(2.2) 
(2.3) 
(2.4) 



Here N 1 is an overlapping factor (|A2|) , L is the distance 
between the mirrors of the FP cavities in arms. We also 
omit non- resonance terms (~ Zi). In order to clarify the 
dependence of zf (A — CI) we have to write down the last 
term in l|2.1|l in detail and to equate exponential quanti- 
ties in the left and right parts of l|2.1() : 
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For the FP cavity 2 all formulas are the same. For the 
mean amplitudes we have: E\ — T\ , £2 = J- 2- 



C Power Recycling Mirror and symmetric mode 
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B. Beam splitter 

We consider that F3, E3 are the field amplitudes on 
the PR mirror and the beam splitter, and F4, E4 are the 
field amplitudes on the beam splitter as it is shown in 
fig. n We assume that the beam splitter transparency 
is Tb s = 1/2 and the phase of the wave due to traveling 
between the FP cavity 2 and the beam splitter is e 1 * 2 = 
1 , and between the FP cavity 1 and the beam splitter is 
e i<t>i _ |_ g Q we can i ma g me such plane (see Fig.QJ, that 
phase advance between the beam splitter and this plane 
is e 1 * 1 = i and between this plane and the input mirror 
of the FP cavity 1 is fold to 2n. Then Fi , Ei are the field 
amplitudes on the input mirror of the FP cavity 1 (and 
on the plane (i)), F2, E2 are the field amplitudes on the 
input mirror of the FP cavity 2 and on the beam splitter. 
For the mean amplitudes we have: 

= -F3/V2, T 2 = -F3/V2, £4 =0, £3 = ^3. 

It is convenient to introduce symmetric and anti- 
symmetric modes Y±: 

F+ = (Ft + F 2 )/V2, T+ = -T 3 , f+ = -f 3 , (2.5) 
F = (Fi - F 2 )/V2, T- = 0, f_ = -if 4 . (2.6) 

Then for fluctuation fields e 3 and 64 we have: 

e 4 =-K(a.)U-Fn n T n V2Nikz*_, (2.7) 
e 3 =n(n)U+^n n T n V2iN^z* + , (2.8) 
z- = (xi -yi) - [X2 -V2), z+ = (xi -yi ) + (x 2 -yi) 

We see that symmetric and anti-symmetric modes can be 
analyzed separately. In subsection III CI we consider the 
symmetric mode which interacts with sum coordinate z + 
and in subsection III Dl the anti-symmetric mode inter- 
acting with differential coordinate z_ . It is worth to note 
that such consideration is possible only if all four mirrors 
in the FP cavities in arms are optically and elastically 
identical, just as we have assumed. For the opposite case, 
when the eigenfrequencies of elastic modes are different 
we can consider elastic mode only in one mirror and as- 
sume other mirrors to be fixed (see subsection III Ef) . 



C. Power Recycling Mirror and symmetric mode 

Assuming that PR cavity is in resonance, we obtain 
two equations in time domain describing coupling be- 
tween optical and elastic modes (see details of calcula- 
tions in Appendix [BJ) : 
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Here yo+ is the relaxation rate of symmetric mode (|B6ll . 
y m is the relaxation rate of elastic mode, m is the mass of 
each of the mirrors, and \± is normalizing constant l)A10(l . 

We find the solutions of this set of equations (|2.91 12.10(1 
in the form: 
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and as a result we obtain the characteristic equation: 
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(2.12) 
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(2.13) 

Here Ai is general overlapping factor between the main, 
Stokes optical modes, and elastic mode, uj_ is the compo- 
nent of displacement vector u of elastic mode normal to 
the mirror surface, J dfj_ corresponds to the integration 
over the mirror surface and J df: to the integration over 
the mirror volume V. 

We present A as a sum of real and imaginary parts: A = 
a+ib, substitute it into l|2.11|l and obtain two equations: 

a 2 -b 2 + a(y m +yo+) + bA + y m y + - 2Q = , 

(2.14) 

2ab + b(y m +Yo+) -aA-y m A = 0. (2.15) 

The parametric instability condition corresponds to a > 
0. Putting a = we find the condition of boundary situ- 
ation (between stability and instability). The additional 
analysis gives the sign of inequality in the parametric in- 
stability condition. We find formal solution for b from 
(|2.15|l with assumption a = 0, substitute it into H2.14f> . 
and obtain the parametric instability condition: 
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This condition can be compared with the condition of 
parametric instability for single FP cavity with relax- 
ation rate y 
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We see that conditions (|2.16|l and l|2.17|l approximately 
coincide with each other if one substitutes yo+ instead 
of y and takes into account inequality Ym *C Yo+ (see 
estimates in Appendix |EJ). The factor 2 appears due to 
the fact that in derivation (|2.16fl we take into account the 
displacements of 4 mirrors in the interferometer while in 
(I2.17|l the displacement of only one mirror of FP cavity 
is taken into account. 



C Power Recycling Mirror and symmetric mode 
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D. Signal recycling mirror and anti-symmetric 
mode 



We assume that in general case the SR cavity is not 
in resonance and phase advance 4> of the wave between 
the SR mirror and the beam splitter is an arbitrary one. 
We also assume that (J> does not depend on frequency 
Q. due to the short length of the SR cavity. After the 
calculations presented in Appendix IClone can obtain two 
equations describing the coupling between optical and 
clastic modes: 



fin- + (TO- — i5) fin- 



NlfHjU,,Z -e^, (2.18) 



LV2 



z_ + y m z_ 



7Ticu m mu 



(2.19) 



Here yo— is the relaxation rate of anti-symmetric mode 
(|C6J| . 6 is the detuning depending on the SR mirror po- 
sition (|C5)l . 

Searching for the solution of the equations set l|2.18l 
I2.19|) in the form: 



fi 



(t) = fi 



(A-iA)t 



t )=z*_e 



A I 



we obtain the characteristic equation: 

2Q= (A+yo--i(6 + A))(A + y m ). 



(2.20) 



We see that the characteristic equation (|2.2U|I differs 
from the analogous one for the symmetric mode (|2.11|) 
only by replacement of A — > A + 6. Hence we can write 
down the condition of instability for the anti-symmetric 
mode using l|2.1ul) : 



2Q 



YmYO- 



> 1 + 



(A + 6) 



(Ym+YO-)' 



(2.21) 



The relaxation rate Yo- and detuning 5 depend on 
angle 4> (the position of the SR mirror). Analyzing def- 
initions (|C6I IC5(I we see that there are two cases: close 
to resonance (4> close to 0) and to anti-resonance (ct> is 
close to 7t/2) cases. Expanding the denominator in (|C6I 
IC5(1 in series over T ST (recall T ST <C 1), we obtain useful 
formulas: 



y sin 2cf) 
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(2.22) 
(2.23) 



Manipulating the angle it) (via the SR mirror position 
variation) we can vary the relaxation rate yo- and de- 
tuning 6 in wide ranges: 



TsrY . 

— < Yo- 



< 



4y — 2y 



' ST 



sr 



< 6 < 



2y 



sr 



(2.24) 




Jo-/y 

5/y 

FIG. 2: Dependence of the relaxation rate Yo- and detuning 
6 of anti-symmetric mode on angle cb for T ST = 0.07 (planned 
in Advanced LIGO). Top: close to resonance case. Bottom: 
close to anti-resonance case. 



It is demonstrated in Fig. [21 (Note that the above for- 
mula is valid only while yo- <C c/L, i.e. T/8T sr <C 1.) 
Using formula (|2.24|) . wc have the following estimates for 
parameters of Advanced LIGO (Appendix |E|) : 

1 .6 sec -1 < y - < 6 x 1 3 sec -1 , 
-3 x 10 3 sec -1 < 6 < 3 x 10 3 sec -1 . 

So we can "scan" the frequency range to find instability 
(or its precursors) by variation of the SR mirror posi- 
tion. As precursors we may register Stokes modes provid- 
ing information about the resonance frequencies of elas- 
tic modes which "suit" each other by spatial distribu- 
tions. These Stokes modes may be the modes of higher 
orders (dipole, quadrupole and so on). In order to ex- 
tract them one has to detune the output mode cleaner 
which is planned to be placed after the SR mirror which 
is an additional but not complicated operation 1 . It pro- 
vides us in situ with very valuable information about the 
possible danger of parametric instability. 



1 Advanced LIGO currently plans a rigid short (tens of cm) out- 
put mode cleaner which can be put on resonance for the car- 
rier but which rejects all frequencies out of mode cleaner band- 
width. Detuning, for example, to dipole mode requires change of 
mode cleaner length i by relatively small value bl/i = (a>di po le — 



E Different mirrors: only one mirror is in resonance 
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E. Different mirrors: only one mirror is in 
resonance 

Here we consider the case when the frequencies of elas- 
tic modes in different mirrors do not coincide with each 
other and thus, we assume that the frequency of elastic 
mode of only one mirror is in resonance, for example, the 
end mirror in the FP cavity 1. Therefore we consider 
other mirrors as fixed, and only the coordinate X] will be 
taken into account. Then Eqs. (|2.91 12.18(1 will be valid 
with substitution X] — > z + , z_. For calculation of pon- 
deromotive force we consider the field in the FP cavity 1 
as a sum of symmetric and anti-symmetric modes fields. 
After calculations presented in Appendix [pj we get the 
characteristic equation in the compact form: 



(A + iA + y T 



--(— ~ 
2 V*+Yoh 



+ 



1 



A — iS +yo- 



(2.25) 



Instability condition in pure power recycled LIGO in- 
terferometer. Using characteristic equation (|2.25|) we 
consider the case when there is no SR mirror in the in- 
terferometer, then we have to substitute T ST = 1 into 
(|C5IC6|) and hence 6 = 0, jo- — J (y is the relaxation 
rate of a single FP cavity in arm). For such power recy- 
cled configuration we have the following inequalities: 



Ym < Yo+ < Y 



(2.26) 



The analysis presented in Appendix[D]gives the following 
parametric instability condition: 



Q ( Yo+ Y 



2ym \Yo + + A 2 Y 2 + A 2 



> 1 



(2.27) 



It is useful to rewrite approximation of this instability 
condition in particular cases for different detunings A: 



A < Yo+ 
Yo+ < A<y 
A>y 



Q 



2jmY0 

Q / Yo+Y 
2y m y V A 2 

Qy 
2y m A 2 



>1, (2.28a) 
+ l) > 1 , (2.28b) 
(2.28c) 



Note that condition (I2.27|) slightly differs from the for- 
mula (6) in our paper Q 2 - However, the error made in Q 
is insignificant: the particular cases l|2.28|l coincide with 
analogous particular cases obtained from formula (6) in 



2 Our notations relate to the notations in as follows: y => 
6l i To+ =¥ 5pr , Tm =r> 6m, and parameters Q and IZo relate as 



The instability condition for signal recycled LIGO in- 
terferometer. Using estimates of Appendix [E] we have 
the inequality 



Ym < Yo+ < Yo- • 



(2.29) 



From the characteristic equation l(2.25|l and inequality 
we obtain the parametric instability condition: 



YCH 



+ 



Yo- 



2y m VY 2 ^ + A 2 y 2 _ + (A + 6) 2 
See details of calculations in Appendix IP 



> 1 . (2.30) 



III. DISCUSSION 



For our discussion wc use the following scale of re- 
laxations rates (see Appendix [EJ): the relaxation rate of 



.-i 



the 



elastic mode Y m ~ 6 x (10~ 4 -r 10~ 2 ) sec 
laxation rates of symmetric and anti-symmetric modes 
Yo+ — 1.5 sec -1 , Yo- > 2 sec -1 , and the relaxation rate 
of a single FP cavity in arm y ~ lOOsec -1 . 

First of all we see that the main difference between 
the parametric instability in a signal recycled Advanced 
LIGO interferometer and in pure power recycled inter- 
ferometer is the crucial dependence on detuning, com- 
pare Eqs. (|2. 1612. 2112. 3U[) with (|2~27|) . In power recycled 
interferometer the parametric instability takes place if 
|A| < y ~ 100 sec -1 while in signal recycled interferome- 
ter: |A| < yo+, Yo- — 2 sec -1 . 

On the one hand, in case of relatively small detuning 
the parametric instability in a signal recycled interferom- 
eter takes place at relatively low value of optical power. 
For example, if A <C Yo+ and 6 ^> yo- one can obtain 
from Eq. (|2.3U|) that the parametric instability will take 
place at power W c ~ 5 W (!) circulating in arms (if 
cu m = 10 5 sec -1 , y m = 6 x 10~ 4 sec -1 , Ai ~ 1). On the 
other hand, there is a small chance that such small de- 
tuning will take place and for large detuning (|A| > Yo+) 
the realization of parametric instability requires dramat- 
ically larger optical power: W c ~ A 2 /yq + . For example, 
if detuning is about 1 kHz (|A| ~ 6 x 10 3 sec -1 ) and other 
parameters are the same one can obtain W c ~ 10 s W (!). 
For the same reason the possibility that the presence of 
anti-Stokes mode can depress the parametric instability 
is small enough especially for such detunings. Therefore, 
we did not consider the anti-Stokes mode in our analysis 
("anti-Stokes" generalization can be done using the same 
approach as in Q). 

Another factor that can decrease the possibility of 
parametric instability is the small value of overlapping 
factor; even if the frequencies of Stokes and elastic modes 
"suit" each other their spatial distributions (at mirrors 
surface) may considerably decrease the overlapping factor 
Ai and, hence, the possibility of parametric instability. 

The elastic modes can be calculated numerically. Un- 
fortunately, the numerical calculations of elastic modes 
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of cylinder mirrors face with obvious difficulty: the accu- 
racy of elastic mode frequencies calculations is dramat- 
icaly insufficient. The standard pac kages FEMLAB or 
ANSYS, used for this purpose [§. Ill|. provide the accu- 
racy about several percents only, while we need the accu- 
racy at least yo+/wm — 10~ 7 10~ 5 (!). Nevertheless, 
numerical calculations have sense for estimates of over- 
lapping factors and as information (about the frequency 
range) for the experimentalist where the parametric in- 
stability may occur. 

Even if the numerical calculation methods will be im- 
proved to the extent that they will achieve the required 
accuracy, they can not solve the problem completely be- 
cause: (i) the shape of mirrors differs from the cylinder 
shape, for example, the pins to attach suspension fiber 
may produce the shift of clastic mode frequency up to 
100 sec -1 Q and (ii) the inhomogeneity of fused silica 
(the material mirrors should be manufactured of) may 
also provide an uncontrollable shift of elastic mode fre- 
quency [U. That is why we present the separate con- 
sideration of the case when the elastic mode of only one 
mirror is taken into account in subsection III El 

Parametric instability can also be investigated in 
GEO600 configuration interferometer which also has sig- 
nal and power recycling mirrors, but has no FP cavities 
in arms. Hence the formulas above can be generalized 
for GEO600 configuration by putting the transmittancc 
of input mirrors in arms to be equal to 1 and, hence, 
y = c/2L. Then the relaxation rate of symmetric mode 
will be about yGEO _ T pT c/4L ~ 0.75 x 10~ 5 sec" 1 (we 
assume T pT = 0.012, L = 1.2 km) and using formula 
(|2.16() we estimate that parametric instability may take 
place at relatively small optical power W^ EO ~ 100 W 
circulating in arms if we assume zero detuning A = 0, 
m = 10 kg, cu m = 10 4 sec -1 , y m = 10~ 4 sec -1 and 
At =1. 

Our consideration can be generalized for the case of 
mirrors with losses. Let each mirrors in the FP cavity to 
have loss coefficient A which is small compared to trans- 
parency: T) = 2A/T < 1. In this case all formulas i|2.16l 
12.211 12.271 I2.3U|) for parametric instability conditions are 
valid with the following substitutions: 

Y=^Y(1+r|), To±^To±+riy. (3.1) 

For Advanced LIGO the losses in mirrors are small 
enough (it is planned that A ~ 5 ppm, n ~ 2 x 10 ). 

However, the diffractional losses may be large for op- 
tical modes with high indices @ . Kip Thorne 0, 0] 
pointed out that the case of large diffractional losses 
(when round trip relative losses are close to unity) re- 
quires a separate analysis. This work is in progress now 



Conclusion 

We have shown that in the signal recycled Advanced 
LIGO interferometer the possibility of falling into the 



trap of the parametric instability is smaller than for the 
pure power recycled one due to stronger sensitivity to 
detuning. 

We think that the most reliable method to avoid the 
parametric instability is the direct experimental test. For 
signal recycled interferometer we have good method to 
investigate the possibility of parametric instability ex- 
perimentally in situ: varying the SR mirror position, 
one can detune the frequency of anti-symmetric mode 
in wide range to find the instability or its precursors as 
it was shown in subsection III Dl It is important that 
for GEO600 configuration we can introduce detuning 
larger than in Advanced LIGO: the scanning may be per- 
formed inside the free spectral range c/2L. This scanning 
combined with the detailed knowledge about the elastic 
modes (it can be obtained in situ in separate experiments 
before the test masses are placed into the interferometer) 
will provide us with very valuable information helpful for 
avoiding the parametric instability. 

We hope that parametric instability effect can be elim- 
inated in Advance LIGO interferometer after the de- 
tailed experimental investigations supported by theoret- 
ical analysis. 
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APPENDIX A: FP CAVITY WITH TWO 
MOVABLE MIRRORS 

Here we derive formulas (12.11 12. 2|) and also obtain the 
equations describing opto-elastic coupling in FP cavity 
with two movable mirrors (see notations in fig. |3J) . 

We denote distance between FP cavities mirrors as L. 
Close to resonance case it is convenient to introduce "gen- 
eralized transparency" 7q and "generalized reflectivity" 
TZn for Stokes mode l|2.3|l using the following approxi- 
mation: 

= e i(tu ' +a)T ~1+iOT, t=-. 

c 

Close to resonance we use the following approximation 
Qx < 1 , T« 1. 

The mean amplitude of main mode falling on end mir- 
ror is fiin and we assume it as a constant. 
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FIG. 3: Scheme of FP cavity and notations. Both mirrors can 
move as free masses. 



In general case for fields on end mirror we have 



A 



1^ (nj 



tcu, t _ _ y~ 



•"Min An) _-iaii t 



Ah 



Vs 



^"une" 1 - a "' t 2ikuj. (xe~ 



Here sum is taken over the complete set AV^ of cavity's 
modes (they are orthogonal to each other) and u_l is 
normal to surface component of displacement vector u 
of elastic mode and x is slow amplitude of displacement. 
Multiplying this equation by distribution function A^ irL 
of our Stokes mode, integrating over cross section and 
omitting non-resonance term (~ xe~ XWmt ) one can find 
in frequency domain 

e U n(n) = -f ni x(0)-N 1 JP- 1in 2ikx*(A-Q), (Al) 
J A 0in A* Un u ± dfi 



N, = 



vSoSi 



(A2) 



For fields on input mirror we have (see notations in 

fig. El) 

eun=iTfi -Rfun + Rt-^unjN^iky*, (A3) 

ei =iTf Un -Rfi -RN^unliaJivVc, (A4) 

fiin = e UT1 .e, f lin = 9eun- (A5) 

Combining these equations one can obtain: 



flin = fl 



iT 



2ik(x*-u*), (A6) 



flin. — Sun.) flin — Clin — — &Mn, (A7) 

ei =f 1 ^a-N 1 ^ Hn T a 2ik(x*-y*). (A8) 

One can find equation for amplitude fun iu time do- 
main using the rule (— iCl) — > 9l: 



f 



lin 



yfun.= ~ 2tk(x*(t) -y*(t))e 



(A9) 



The light pressures acting on end mirror and input 
mirror inside the cavity are approximately equal to each 
other. So we can calculate light pressure acting only 



on end mirror. Keeping only cross term we can obtain 
formula for light pressure P in time domain: 



-^oin^Hn-^Hnfinl M e 



-i(o>o-o>, )_ 



cvSoSi 

4- A* A,- T* f. -, f+1 p' 1 ( cu o-"Ji ) 



Now we can write down equation for elastic oscillations 
with amplitude x: 

p Y_ u« (x* + 2y#xt + (o># ) 2 x £ ) = rt ± V(f ± ) 6(T|| - rjj) . 



Here p is density of mirror, sum is taken over the com- 
plete set of elastic modes (spatial displacement vectors 
U£ are orthogonal to each other), n± is unit normal to 
mirror's surface, tm is longitudinal coordinate of points 
inside body of mirror, coordinate rjj corresponds to face 
surface of mirror. Multiplying this equation by spatial 
distribution vector u* of our elastic mode and integrat- 
ing over mirror volume V one can obtain: 



x + 2y m x + u> m x = 



2(N?J7 ta fim(t)e 



iltuo-tu, H . 



emu. 



u(f)| 2 df, 



(A10) 



where m = pV is mirror's mass. 

Introducing the slow amplitudes x xe~ tcUmt + 
x*e lWmt we write down equation for amplitude x*: 



N 



tcmucum 



J7 tn firv1 (t) e 



iAt 



(All) 



And for the coordinate z = x — y we finally obtain (if FP 
mirrors are elastically identical): 



2N 



YmZ 



tcmutum 



^ 4n fmi (t) e 



iAt 



(A12) 



APPENDIX B: SYMMETRIC MODE 

In this Appendix we derive equations l|2.9l I2.10fl for 

analysis of parametric instability in symmetric mode. 
One can start from equations for amplitudes F3 and E3 
on beam splitter: 



VWs- V^Ese 1 *", (Bl) 
i/VEse 1 *- -V^~ V F 5, (B2) 

(cu + Apr +Q)lp r /C. (B3) 



We assume that PR cavity is in resonance: exp(icb pr ) = i 
and we assume that (\> vr does not depend on frequency 
CI due to shortness of PR cavity (l pr <C L). Then using 
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(12.8(1 one can obtain: 

y — in — f 5 



£5 



(yo+ - iQ) 1 + ^1 — Tp T 

A/T(l + /^V)( Yo+ -in)' 
f 5 (Yo++in) \/T pT ^unY2V2N 1 kz* 



(B4) 



(Yo+-iO) VT(1 + ^1 -T pT )(Yo+-ifl)' 

(B5) 



Yo+ =Y 



1 V^l Tpr ^ TprY 



(B6) 



XTie fluctuation part f + o/ symmetric mode. We 
rewrite equation l|2.1|l using (|2.5IIB4|) : 



fin+ — 



2iVTp r Y f 5 



VT(1 V)(Yo + -iO) 
^ Un 2x/2iYlcNiz; 



+ 



(B7) 



+ 



T(Yo+-in) 



Recall that in frequency domain we mean z* + — z* + (A— CI) 
(A is a detuning (E2J)- 

Now we can write down Eq. 1B7|) in time domain using 
rule — iCl — > 3 t : 



i-v/VV f 5 , 2V2-uFi in N 1 Ykz* 
(3t +Yo+Jfin+ = — — t= + 



VT 



T 



For analysis of parametric instability in last equation we 
omit term proportional to fs and finally obtain Eq. I|2.9[) . 

Then we can calculate the pondcromotivc forces acting 
on each mirror and write equation for sum coordinate z + : 



pm- 

m 



1 pm.4 



z+ +2y m z+ + w m z+, 



2V2Ni 



7t 



,i(a)o— o>i 



— i(a> — u>, )_ 



(B8) 
(B9) 



We rewrite |B9j) keeping only resonance terms: 



7t 

-,io) m t / • * 



Z+ + TmZ 4 



(BIO) 
(BH) 



2icu m m 

Now we can obtain equation l|2.10[l . 

APPENDIX C: ANTI-SYMMETRIC MODE 

In this Appendix we derive equations (|2.18l I2.19|) 

for analysis of parametric instability in anti-symmetric 



mode. For amplitudes F4 and E4 on beam splitter we 
have: 



i/f^Fg - VTV^e 1 *, (CI) 
i/f^E 4 e i * - v/^T^^e. (C2) 



We assume that SR cavity is not in resonance (i.e. 4> = 
(cuo + Cl)l ST /c is an arbitrary number) and that 4> does 
not depend on frequency II due to shortness of SR cavity 
(l sr <C L). Then using (|2.7(l one can obtain: 



U = 



if 6 VW*(Y-ia) 



;i + VT~TTr e 2i *) (yo- - i(5 + CI)) 
| 2v / 2iYe 2ic >/r~T^^unNikzl 

VT(1 + yi~n^e 21 *) (yo- - 1(5 + CI)) ' 
fc e 2i * (y -+i(5 + Q)) 



(C3) 



(C4) 

T(1 + VT~T7 r e 2i ^) (yo- - i(6 + a)) ' 

(C5) 



(yo_-i(6 + Q)) 

2V2 T VT sr c l *^unNilcz* 



Y\/1 — T ST sin 2ct) 



Yo- 



(1 -T^/2 + 71 - T^cos24))' 

Yjsr 

2(1 -T sr /2 + vT^: cos 2cb)' 



(C6) 



We rewrite equation 112. 1[) using l|2.6l IC3|) to obtain 
formula for fluctuation part f_ of anti-symmetric mode: 



2YV /= f^fee i * 



VT(1 + /I - L7e 2 ^) (yo- - i(6 + O)) 



(C7) 



2V2iY^n^Nikzl 
T(y --i(S + Q}) 



For analysis of parametric instability we can omit the 
term proportional to f s so that 



fin- = 



NiF nn iwiz*_ 
lV2(Yo--i(6 + n)) 



(C8) 



Finally one can obtain from (|C8|) the equation 12.1811 in 
time domain. 

In the same manner as for symmetric mode one can 
obtain the formula for ponderomotive force in anti- 
symmetric mode: 



pm— . ~ 

= z_ + 2y m z 

m 

2V2S 



p ra- 



ft 



^linfir.-We 



-1 U)n— O), t 



(C9) 
(CIO) 



^T ln fin-(t)e 



ilujQ-UJ! )t 



Keeping only resonance terms and rewriting (|C9IC10|) 
one can obtain equation l|2.19|l for differential coordinate 
z_ in time domain. 
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APPENDIX D: ONLY ONE MIRROR IS IN 
RESONANCE 



For ponderomotive force we have 



2S(E 



main ~t~ t-Stokcsj 

4n 



(Dl) 



(Nl ^un(t) fJnl (t) + N^ 4n (t) fin! (t)) 



Keeping only resonance term (~ fi n i ) we obtain equation 
for elastic mode: 

N*J-"*- e iAt 

Xi+YmXf^^- 1 — — flin(t), A = U) - U)l - W m 

icnmcu m 

Taking into account both anti-symmetric and symmet- 
ric modes we finally obtain the set of equations in time 
domain: 



*1 + Ym*1 



^T-^fin Hn+ + fin- iAt 



V2 



0, 

(D2) 



Xi e 



-iAt 



(3t +T0+)fin+ =0, 



V2L 



Xi e 



-iAt 



(D3) 

(9 t -i6 +yo-)fin- =0. 

(D4) 

We sum and subtract equations (|D3I ID4|I (see also defi- 
nitions 

r - 



; . Yo+ + T - To 

Min H ^ THn H 



; . Yo+ + T - 

T2in H ~ T2in 



2 



2 2 

where we introduced notations: 



f2ln = ff, (D5) 
f U n=0, (D6) 



iNl.Fl in "J 



1 i-At 



Manipulating Eqs. I|D5I ID6|) we get: 



3t + 



Yoh 



x ID 



Yoh 



r - 



flin + IY0+ + To-Jflin + YO+To- THn = 0" H ^ CT " 

Finally we have the set of equations: 

N* f* 

** + Ym*i - fun(t)e 1At = 0, (D7) 



-iNiFu n o>i / 9 + y 



icnmcu m 

r - 



(D8) 



L ' 2 

+flin + (Y0+ + T - )f Hn +Y0+T - f Hn = . 

Finding the solution of this set in form: fun(t) = 
fiine At , x*(t) = x |e (A+lA)t we obtain characteristic 
equation: 

(A + iA + Ym) (A + Yo+) (A + r _) J2 
(2A + Yo+ + r -) ~2" 
Using this equation we get Eq. i|2.25[l . 



TTie instability condition for pure SR configuration. 
As before we substitute A = a + i(b — A) into (|2.25|) 
taking a = 0: 



Ym =Yo+^ + Yo-A2, 

Tml Ytn2 

b = Ai (A - b) + A 2 (A + 6 - b) . 



Here we introduce notations 
2 1 



At 

A 2 



2y2 + + (b-A) 2 
Q 1 



, Yml =Yo+Ai, 



(D9) 
(D10) 

(Dll) 



, Ym2 = Yo-A 2 . (D12) 



2 y 2 _ + (b - A - b) 2 
We can formally solve Eq. (|D10|) and find b 



I^i A +^(A + 6) 



Yoh 



Yo- 



1 + Tm1 , Tm2 



Yo+ Yo- 



(D13) 



Note that values y m i and y m 2 in (|D9|) are positive ones 
and hence Ymi,Ym2 < Ym- Also taking into account 
that Ymi ,Ym2 -C Yo+,Yo- we can conclude from HD13|) 
that 



b < |A|, |A + 5| 



(D14) 



Then the parametric instability condition l|2.3L)[l can be 
easily obtained from Eq. (|D9|) using inequality l|D14(l . 



APPENDIX E: NUMERICAL PARAMETERS 



x JD6J : We used the parameters planned for Advanced LIGO 



5] 



m 
T 

Y 

T sr 



2 x 10 15 sec" 1 
40 kg, 
5 x 1(T 3 , 
94 sec" 1 , 
7 x ICr 2 . 



L 
W 

Yo+ 



4 x 10 5 cm, 
830 kW, 
6 x lO- 2 , 
1 .5 sec -1 , 



We also assume that FP mirrors are fabricated from 
fused silica with angle of structural losses $ — 1 .2 x 
10~ 8 and for elastic modes frequencies cu m = 10 5 -j- 
1 7 sec -1 we estimate relaxation rate Ym of elastic modes 
by formula: 



Tm = w m 4>/2~ 6 x (10~ 4 -7- 10~ z )sec 



A-2 \„„„-1 
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